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Abstract. We discuss recent developments in assessing heavy-quark interactions in 
1—5 ■ the Quark-Gluon Plasma (QGP). While induced gluon radiation is expected to be 

^\ \ the main energy-loss mechanism for fast-moving quarks, we focus on elastic scattering 

which prevails toward lower energies, evaluating both perturbative (gluon-exchange) 
and nonperturbative (resonance formation) interactions in the QGP. The latter are 
| treated within an effective model for D- and B-meson resonances above T c as motivated 

t— I ■ by current QCD lattice calculations. Pertinent diffusion and drag constants, following 

from a Fokker-Planck equation, are implemented into an expanding fireball model for 
' Au-Au collisions at RHIC using relativistic Langevin simulations. Heavy quarks are 

\Q . hadronized in a combined fragmentation and coalescence framework, and resulting 

electron-decay spectra are compared to recent RHIC data. A reasonable description of 
<i ' both nuclear suppression factors and elliptic flow up to momenta of ~ 5 GeV supports 

the notion of a strongly interacting QGP created at RHIC. Consequences and further 
f*> |' tests of the proposed resonance interactions are discussed. 

<D '. 
43 ■ 

> : 

^ . 1. Introduction: The Virtue of Heavy Quarks in Heavy-Ion Collisions 

The observed suppression of high-momentum particles in central Au-Au collisions is one 
of the most exciting discoveries made at the Relativistic Heavy-Ion Collider (RHIC) to 
date. It has been interpreted as a large energy loss of energetic partons traveling through 
an almost opaque partonic medium, the "strongly-interacting Quark-Gluon Plasma" 
(sQGP). The underlying microscopic mechanism has been attributed to medium-induced 
radiation of gluons. However, recent data at RHIC, showing a strong nuclear suppression 
(small Raa) [H 12] and elliptic flow (large v^) [3] of "non-photonic" single electrons 
(attributed to the semileptonic decay of charm and beauty hadrons) in semi-/central 
Au-Au, have posed challenges to the radiative energy-loss picture within perturbative 
QCD (pQCD). The importance of elastic interactions of heavy quarks (HQs) in the QGP 
in this context has been pointed out in Refs. |H El El 13 IE] • Within pQCD, radiative 
and elastic energy-loss mechanisms turn out to be equally important in the currently 
accessible momentum range at RHIC, and their combined effect reduces the discrepancy 
with the measured electron Raa in central Au-Au (but does not eliminate it) [Zj. 
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However, even when using upscaled transport coefficients within pQCD energy-loss 
calculations the maximal single-electron (e ) Vi is limited to 2-3% in semicentral Au- 
Au, while the experimental values reach up to 10% around pf*=2 GeV [B] . Significantly 
larger theoretical v 2 values can only be obtained if the quarks interact strongly enough 
to become part of the collectively expanding medium, implying that not only energy- 
loss but also energy-gain processes (detailed balance) need to be accounted for. This 
has recently been studied employing Langevin simulations of HQs in an expanding QGP 
based on elastic reinteractions jHEJEEIBj- Since at RHIC energies, secondary production 
of HQs is suppressed jT0"] . the latter are also valuable probes of the soft (collective) 
properties of the putative sQGP, and can provide tests of hadronization mechanisms, 
e.g., quark coalescence [TI]. HQ-momentum distributions furthermore have important 
impact on secondary production (regeneration) of heavy quarkonium states [T2 | IT3 t IT^j . 

In the present paper, we report on our studies |U Ej of HQ properties in the sQGP 
and applications to RHIC, based on non-perturbative interactions. We will address 
pertinent uncertainties, open problems and consequences of our findings thus far. 

2. Heavy-Flavor Baseline Spectra at RHIC 

The interpretation of semileptonic e ± -spectra in A- A collisions requires a reliable 
estimate of baseline spectra in p-p, in particular their decomposition into contributions 
of charm and bottom decays (bottom quarks are much more inert to changes in their 
momentum spectra). Current pQCD predictions expect the transition from charm to 
bottom at around p e T ~4-5GeV, albeit with significant uncertainty jTH]. This rather 
early onset of a bottom component constitutes part of the difficulty in explaining the 
large suppression observed in the e ± -spectra in Au-Au collisions. In Refs. jHl QSj, we 
have adopted a model-independent procedure to determine e ± baseline spectra, by first 
reproducing available data on D and D* spectra, evaluating their contribution to the 
e ± spectrum at low(er) p?, and then adjusting the bottom contribution to fill the high- 
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Figure 1. Left panel: D and D* p^-spectra in v / s]vjv=200 GeV d-Au collisions |17|: 
right panel: composition of semileptonic electron-decay spectra based on the charm 
spectra in the left panel in p-p and d-Au at RHIC[TS). 
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Figure 2. Left panel: mesonic spectral functions in the qq channel for "strange" 
quarks note the approximate degeneracy of the chiral multiplets. Right panel: 
total cross sections for a c quark scattering off partons within LO QCD (blue curves), 
and off antiquarks in the resonance model (red curves). 



Pt part of the spectrum, cf. Fig. ^ Our result confirms the crossing of the c and b 
contributions at about pr— 5 GeV. 

3. Elastic Heavy-Quark Scattering in the QGP 

We evaluate (elastic) interactions of HQs in the QGP using Brownian motion of a 
heavy particle in a thermal background of light partons, as encoded in a Fokker-Planck 
equation [13 1201 HUH, 

dt dp dp 2 

for the HQ distribution function, /. 7=Tq 1 and D p are the corresponding drag and 
(momentum) diffusion constants which determine the approach to equilibrium and 
satisfy the Einstein relation, T = D p /^Mq. They are typically calculated from 
elastic 2^2 scattering processes, p + Q — > p + Q (p=q, q, g; Q=c, b). In leading- 
order (LO) pQCD these are dominated by t-channel gluon exchange (regulated by a 
Debye mass). E.g., at a temperature T=300 MeV, and for a strong coupling constant 
a s =0A, the thermal relaxation time for charm is r c ~15 fm/c, well above expected 
QGP lifetimes of tqqp<5 fm/c at RHIC. On the other hand, lattice QCD suggests 
that hadronic resonance (or bound) states in both Q-Q and q-q channels might survive 
up to temperatures of ~2T C j^U 122] , cf- left panel of Fig. |21 In Ref . jl] we therefore 
suggested that D- and I?-meson resonances in the Q-q channel could accelerate the 
thermal relaxation of HQs. Based on an effective Lagrangian for the q-Q-<& interaction 
($=D, 5), 

£ = Q^$rg + h.c (2) 

we evaluated elastic Q + q — >• Q + q scattering amplitudes via $ exchange in s-, t and 
w-channel. Assuming the existence of one $ state (e.g., a pseudoscalar J p =0~), and 
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Figure 3. Left panel: thermalization times of c quarks in a QGP using LO pQCD 
(upper band) and resonance interactions (lower band); right panel: c- and 6-quark 
thermalization times for LO pQCD interactions (a s =0.4, upper green and blue curve) 
and when adding resonance interactions (r<j,=0.5 GeV, lower green and red curve). 



a minimal degeneracy following from chiral and HQ symmetries (with Dirac matrices 
T=l, 75, 7 M , 757 M in Eq. (j2j)), the tq's are reduced by a factor ~3 over pQCD scattering 
(cf. Fig. EI), with moderate sensitivity to the unknown coupling at the effective c-q-D 
(b-q-B) vertex. While the total cross sections for LO pQCD and resonance scattering 
are not very different in magnitude (right panel in Fig. EJ), the angular distributions 
are: forward-dominated pQCD scattering is much less efficient in isotropizing the 
momentum distributions compared to the isotropic resonance scattering, since the 
(angular-averaged) "transport cross section" carries a weight of 1 — cos# (6: scattering 
angle in the center of mass). 
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Figure 4. HQ nuclear suppression factor (left panel) and elliptic flow (right panel) 
following from relativistic Langevin simulations for an expanding QGP in Au-Au at 
RHIC; red (green) lines: c (b) quarks using resonance+pQCD interactions; blue lines: 
c quarks with LO pQCD scattering only. 
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Figure 5. Comparison of the c-quark Raa and V2 following from Langevin simulations 
at RHIC when using effective q-c-D vertices with renormalization (green curves) or 
monopole form factors with cutoff A=l GeV (red curves). 



4. Heavy-Quark and Single-Electron Spectra at RHIC 

Using the initial c- and 6-quark spectra as discussed in Sec. the diffusion and drag 
coefficients as evaluated in Sec. E] have been implemented into a relativistic Langevin 
simulation for an expanding QGP fireball at RHIC 6 j. As expected, resonance 
interactions have a much larger effect on the nuclear suppression factor (Raa) and 
elliptic flow (v 2 ) of c quarks in semi-/central Au-Au collisions than pQCD interactions 
(cf. Fig. EJ), while b quarks are less affected. The leveling-off of the c-quark elliptic flow 
characterizes the transition from a quasi-thermal to a kinetic regime. 

Fig. 03 illustrates the sensitivity of the Langevin results to the regularization 
procedure applied to the Q-q loops [4|: at quark momenta above 3 GeV a renormalized 
interaction (point-like vertices) leads to somewhat stronger effects compared to finite- 
size vertices (form factor with cutoff A=l GeV and identical resonance widths). 

To compare to observables, the HQs have been hadronized and decayed 
semileptonically fTJ E]- If the hadronization is performed by using 5-function 
fragmentation only, the resulting e ± Raa and v 2 do not describe the data very well, 
cf. Fig. El The situation improves when allowing for coalescence of the HQs with light 
quarks fH\ E] at T c which preferentially occurs at lower momenta (the remaining HQs 
are again 5-function fragmented). Consequently, we find an increase in both Raa and v 2 
mostly for momenta up to p^~4-5 GeV, leading to fair agreement with experiment, 
cf. FigEl At higher p?, contributions from radiative energy-loss (which have been 
neglected here) are expected to become important. Note that at high pt the data 
for Raa stay small, while the v 2 tends to (much) reduced values, which is suggestive 
for substantial energy loss with small collectivity. Such a feature is indeed consistent 
with perturbative energy-loss mechanisms. We also recall that the Langevin results in 
Figs. El and presumably overestimate the suppression and v 2 at higher somewhat, due 
to the assumptions of (i) point-like vertices (cf. Fig. EJ), and (ii) a thermal background 
medium even at high p? (while from light hadron spectra of RHIC one infers that full 
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Figure 6. Nuclear suppression factor (left panel) and elliptic flow (right panel) 
for non-photonic single-electron spectra in semi-/central Au-Au collisions at RHIC. 
Data El are compared to theory predictions 6 using Langevin simulations with 
elastic c- and 6-quark interactions in an expanding QGP fireball and ^-fragmentation 
at T c ; red band: LO pQCD + resonance interactions (with r$=0.4-0.75 GeV), blue 
line: LO-pQCD only, purple line: c quarks only for LO pQCD and resonances. 
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Figure 7. Same as Fig. [Tj] but including heavy-light quark-coalescence at 
hadronization. 



thermalization for light partons reaches to about pr— 1 GeV). 
5. Discussion and Further Tests 

The small Raa and the large t>2 i n the e^ data require strong interactions of heavy 
quarks in the QGP, (well) beyond radiative energy-loss. While elastic pQCD scattering 
reduces the discrepancy in the Raa, a sufficient collectivity to explain the t>2 can only 
be obtained for (unrealistically!?) large coupling constants jS]. Resonance scattering is 
a promising candidate for a microscopic explanation of the apparently large diffusion 
constant (for an estimate of energy loss due to 3^3 scattering, cf. Ref. |2~3]). 

An important feature that distinguishes resonance contributions from pQCD relates 
to the chemical composition of the QGP. While the resonances as implemented here 
require the presence of anti-/quarks (here we have assumed QGP in chemical equilibrium 
for 2+1 flavors, amounting to ~25 anti-/quark and 16 gluon degrees of freedom), 
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pQCD approaches typically assume a gluon plasma by saturating the bound on the 
produced entropy in central Au-Au at RHIC with diV g /d?/=i000; this maximizes the 
color charge and thus the interaction strength in both radiative (gluon emission) and 
elastic (gluon exchange) channels. Thus, if RHIC indeed produces a gluon plasma, 
resonance scattering could be severely suppressed (much less is known about the presence 
of colored composites [24]). However, if quark-pair production occurs early in the 
collision resulting in an approximately chemically equilibrated QGP more than 
half of the partons are anti-/quarks and the pQCD transport coefficient is appreciably 
reduced, increasing the importance of non-perturbative effects. 

While the existence of resonances needs to be scrutinized theoretically (e.g., 
by 1QCD calculations of heavy-light correlators, or effective models based on 1QCD 
potentials OEHj), one should also search for experimental means of discrimination. The 
relative quark content of the colliding nuclear system can be increased by lowering the 
center-of-mass energy, while antiquarks (and possibly gluons) become suppressed. Thus, 
one would expect stronger rescattering of c quarks compared to c quarks, resulting in 
larger t> 2 and smaller Raa fc> r D~ and D° relative to D + and D°-mesons. Whether such 
an effect is detectable at lower RHIC energies, or whether one has to await the (higher 
luminosities at the fixed-target) Compressed Baryonic Matter (CBM) experiment at 
the future GSI facility, needs to be checked in quantitative studies. We also note that 
diquark-like c-q interactions could diminish the differences between D and D mesons 
(or even charm-exchange reactions in the hadronic phase |27j). 

Finally, it is important to keep in mind the consistency of open and hidden charm. 
As has been emphasized by several authors |2EJ EE21 EH El, secondary production 
(regeneration) of charmonia is much facilitated by a softening (thermalization) of the 
c- and c-quark momentum distributions (for bottom/onia, this interplay appears to be 
less pronounced, even at the LHC [29J). 

6. Conclusions 

Heavy-quark probes at RHIC have posed new challenges to the theoretical understanding 
of parton-energy loss at high momenta and (the approach to) equilibration at moderate 
and low momenta. Currently, the severeness of this problem hinges on the interpretation 
of (semileptonic-decay) electron spectra, and in particular their relative decomposition 
into charm- and bottom-decay contributions (obviously, explicit D-meson measurements 
would remove this ambiguity). If their crossing in baseline p-p spectra occurs around 
p T ~5 GeV, the combination of radiative and elastic pQCD interactions of heavy quarks 
is apparently too weak to generate a suppression and collectivity consistent with the 
observed spectra in semi-/central Au-Au collisions. We have presented a scenario based 
on resonance exchange to introduce nonperturbative elastic interactions of heavy quarks 
in the sQGP at moderate temperatures (relevant for RHIC). When treated within a 
Brownian motion scheme for an expanding QGP fireball, a fair description of the single- 
e ± spectra up to px— 4-5 GeV emerges (especially if hadronization is supplemented with 
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heavy-light-quark coalescence); at higher p T , induced gluon radiation (not included 
here) is expected to take over and become the main source of heavy-quark energy loss. 
A consistent implementation of the latter in the Fokker-Planck approach, which requires 
the integration of coherence effects, is not an easy task. We have also emphasized the 
importance of the chemical composition of the (early) QGP which affects the diffusion 
constants in the perturbative and nonperturbative sectors (as evaluated here) in a rather 
opposite way and may thus help to disentangle the two. An excitation function could 
provide a valuable handle to explicitly vary the quark-to-gluon ratio to further illuminate 
the nature of heavy-quark diffusion in the QGP. 

Acknowledgment 

One of us (RR) would like to thank the conference organizers for the invitation to a 
very informative meeting. This work was supported in part by a U.S. National Science 
Foundation CAREER award under grant no. PHY-0449489. 



References 



Adler S S et al. [PHENIX Collaboration] 2005, Phys. Rev. Lett. 94 082301 

Bielcik J et al. [STAR Collaboration] 2005 arXiv: |nucl-ex/ 0511005 

Adler S S et al. [PHENIX Collaboration] 2005, Phys. Rev. C 72 024901 

van Hccs H and Rapp R 2005 Phys. Rev. C 71 034907 

Moore G D and Teaney D 2005 Phys. Rev. C 71 064904 

van Hees H, Greco V and Rapp R 2006 Phys. Rev. C 73 034913 



Wicks S, Djordjevic M, Horowitz C and Gyulassy M 2005 arXiv: |nucl-t h/0512076 
Goissaux P B 2006, these proceedings 



Armesto N, Cacciari M, Dainese A, Salgado C A and Wiedemann U A 2005 arXiv hep-ph/0511257 

Levai P, Muller B and Wang X N 1995 Phys. Rev. C 51 3326 

Greco V, Rapp R and Ko C M 2004 Phys. Lett. B 595 202 

Grandchamp L and Rapp R 2002 Nucl. Phys. A 709 415 

Grandchamp L, Rapp R and Brown G E 2004 Phys. Rev. Lett. 92 212301 

Thews R L and Mangano M 2006 Phys. Rev. C 73 014904 

Cacciari M, Nason P and Vogt R 2005 Phys. Rev. Lett. 95 122001 

Rapp R, Greco V and van Hees H 2005 arXiv hep-ph/0510050 

Adams J et al. [STAR Collaboration] 2005, Phys. Rev. Lett. 94 062301 

Suaide A A P et al. [STAR Collaboration] 2004 J. Phys. G 30 SI 179 

Svetitsky B 1988 Phys. Rev. D 37 2484 

Mustafa M G, Pal D and Srivastava D K 1998 Phys. Rev. C 57 889 
Asakawa M and Hatsuda T 2004 Phys. Rev. Lett. 92 012001 

Karsch F and Laermann E 2003 in Hwa R C and Wang X-N (eds.) Quark gluon plasma 1, and 

arXivj |hep-lat/0 305025 
Liu W and Ko C M 2006 arXiv nucl-th/0603004 
Shuryak E V and Zahed I 2004 Phys. Rev. C 70 054507 
Lappi T 2006 arXiv |hep-ph/ 0606090 
Mannarelli M and Rapp R 2005 Phys. Rev. C 72 064905 
Cassing W 2006 private communication 

Braun-Munzinger P and Stachel J 2000 Phys. Lett. B 490 196 

Grandchamp L, Lumpkins S, Sun D, van Hees H and Rapp R 2005 Phys. Rev. C in press, and 
arXiv hep-ph/0507314 



